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Abstract — A model with flat FRW symmetries and G and A, vari- 
able is considered in such a way that the momentum- energy tensor 
that describes the model is characterized by a bulk viscosity parame- 
ter. For this tensor the conservation principle is taken into account. 
In this paper it is showed how to apply the dimensional method in 
order to solve the outlined equations in a trivial way. 

Keywords — Odes, AD, FRW Cosmologies, variable 
canstant 



I. Introduction. 

Recently several models with FRW metric, where "con- 
stants" G and A are considered as dependent functions on 
time t have been studied. For these models, whose energy- 
momentum tensor describes a perfect fluid, it was demon- 
strated that G (xt", where a represents a certain positive 
constant that depends on the state equation imposed while 
A oc is independent of the state equation (see 0,1^). 
More recently this type of model was generalized by Arbab 
(see [^) who considers a fluid with bulk viscosity (or sec- 
ond viscosity in the nomenclature of Landau (see Q)). The 
role played by the viscosity and the consequent dissipative 
mechanism in cosmology has been studied by many authors 
(see [|).). 

In the models studied by Arbab constants G and A are 
substituted by scalar functions that depend on time t. The 
state equation that governs the bulk viscosity is: ^ oc ^gp''' 
where 7 is a certain indeterminate constant for the time 
being 7 e [0, 1]. 

As we shall see, this problem is already solved, but our 
aim is to solve it through Dimensional Analysis. We mean 
to point out how an adequate use of this technique let us 
find the solution of the outlined equations in a trivial way, 
even pointing out that it is not necessary to impose the 
condition div{Tij) = 0. The paper is organized as follows: 
in the second section the model is showed, expounding the 
equations and showing the ingredients that compose the 
model. The third section is devoted to revise the solution, 
which is reached by means of standard techniques of ODEs 
integration (see in special while in the forth section 

the dimensional technique will be developed in order to 
solve the model. This section is divided in two subsec- 
tions. In the first one, titled "pretty simple method", we 
point out how a naive use of D.A. brings us to find such 



solution in a trivial way. Several cases are also studied here 
by Arbab I. Arbab (see §), while in the other subsection 
a finer dimensional technique is showed. That is why we 
call it "not so simple method^ . This section is based on 
dimensional techniques (groups and symmetries, see Q), 
in order to reduce the number of variables intervening in 
the expounded ODEs. They are so simplified that its inte- 
gration is immediate. We think that the technique showed 
here is so powerful that it shall be proved that imposing 
the condition div(Tij) = is not necessary to impose in 
order to solve the equations. 

II. The model. 

This problem was posed by Arbab (see |^). The equa- 
tions of the model are: 



Ri 



-gijR- A{t)gi. 



„4 



(1) 



and it is imposed that^: 



div{Ti, 







where A{t) represent (stand) the cosmological "constant" . 
The basic ingredients of the model are: 
1. The line element defined by: 



ds' 



cMr + f\t) 



dr^ 



1 - kr 



{de 



2 I ■ 2 
sm 



we only consider here the case fc = 0. 

2. The momentum-energy tensor defined by: 



(2) 



where p is the energy density and p* represents pressure 
[p] = \p*]- The effect of viscosity is seen in: 



p* = p~ 3^H 



(3) 



where: p is the thermostatic pressure, H — {f'/f) and ^ is 
the viscosity coefficient that follows the law: 



(4) 
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^we shall see that this condition it is not necessary to impose it 
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where makes the equation be homogeneous i.e. it is a 
constant with dimensions and where the constant 7 G [0, 1]. 
And p also verifies the next state equation: 



p = u)p 



Lu = const. 



(5) 



where u £ [0, 1] (i.e. it is a pure number) so that the 
momentum-energy tensor verifies the so-called energy con- 
ditions. 

The field equations are: 



,r ^ (T) 



l\2 



87rG{t) , 



- c^A{t) 



iff S^Gjt) 

/2 



deriving (7) and simplifying with (6) it yields 

A'c4 G' 



p' + 3(w -t- l)pH - Qk^p'iH^ + 



P- 







SttG G 

and at the moment we consider this other equation 

div{T,j)^Q ^ p' + 3{p + p*)L ^0 

if we develop the equation (9) we get: 

p' + 3(w + l)pH - 9k^p^H^ = 
III. Non Dimensional Method. 



(6) 
(7) 

(8) 

(9) 

(10) 



In this section we will mainly follow Singh et al work (see 
10|). If we take the equation M) regrouped, we get: 



p' + 3(w + l)pH - 9k^p''H^ 



G' 

Pg 



8ttG 



(11) 



Al 



A2 



if take into account the conservation principle 
p' + 3(w + l)pH - 9k^p^H^ = 



(12) 



then we solve this equation by solving the equation A2 in 
(11), in such a way that the equation to be solved is now: 



G' 

Pg 



SvrG 



= 



(13) 



this equation is tried to be solved like this (see [10]). It 

is defined A = "^^^ where /3 is a numerical constant, 
(hypothesis by Arbab (see [^) as well as by Singh et al 
(see |lOj), condition that as we shall see, it is not neces- 
sary to impose in the solution through D.A.) and from 
the equation (7) the following relationship is obtained: 
SirGp = 3(1 — j3)H^ . Hence if all the equalities are re- 
placed in the equation (|l3| ) it yields: 



H' 



{l-(3)H 



(14) 



which is easily integrated. 

H = Cip^/'^ 



(1-/3) 

we get to the equation (12) with all these results 

p' + 3(w + l)pH - 9k^p^H^ = 
we arrive to the next equation: 



(15) 



p + 3Ci(w + l)p— - 9Gfk^p- 







which has got a particular solution in the case 7 
obtaining: 



(16) 



Pit) 



1 



{aot) 

and obtaining from it: 



/ ao = (3Ci(w + 1) - 9k^Cf) d'^ 



/(i)=C2i(3(t^ + l)-3fc^Ci)(l-7) 

This is the most developed solution reached by Singh et al 
(see [10]) which is slightly different from the one by Arbab 
(see i). 

IV. Dimensional Method. 

We shall explore this section two dimensional methods. 
The first one, probably the simplest one, has the incon- 
venience of having to depend on Einstein criterion(see [Q 
and Barenblatt [^), while the second one is more pow- 
erful and more elaborated. We shall finish showing an 
equation obtained without having to impose the condition 
div{Tij) = 0. 

A. Simple Method. 

The dimensional way followed in this section is probably 
the most basic and simplest one. On one hand we integrate 
independently the equation 



f 

div{Tij) = -^p' + 3{lu + l)pij 







(17) 



not taking into account the term 9k^p~' , since if we cal- 
culate its order of magnitude we verify that is very small 
K, 10~*° following, then, an asymptotic method (or pertur- 
bative) but this must be justified from a physical and/or 
mathematical point of view. If we integrate the equation 
(17) it is obtained the well-known relationship: 



-3(w+l) 



(18) 



from this equation it is obtained one of the dimensional con- 
stants of our problem: Ai^, that has different dimensions 
and physical meaning depending on the state equation im- 
posed i.e. it depends on The other dimensional constant 
considered has been obtained from the state equation (4) 
i.e. ^ = k^p'^, such constant k-y will also have different di- 
mensions depending on the value 7, in such a way that the 
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problem is reduced to the following set of quantities and 
constants 971. 

M = {t, c, A^,kj, a) 

where its respective dimensional equations in regard to a 
base B = {L, M, T, 9} are (the base B of this type of mod- 
els has been calculated in 0): 



[t] = T [c] = LT-i 



1 _ ^7-l^,^l-72i2ri-l y^^J 



where a represents the radiation constant and it will be 
take into account when we consider the thermodynamics 
quantities. 

Having done these considerations our aim is, therefore to 
solve this model through D.A. The Pi-theorem will bring 
us to obtain two tt dimensionless monomials; one of them 
will be the obtained in the case of a perfect fluid ( [2] ) and 
the other monomial will contain information on viscosity, 
showing in this way that this type of models are very gener- 
al, reproducing the results obtained in the case for perfect 
fluids. Since all solutions will depend on these two mono- 
mials we must take into account Barenblatt criterion if we 
mean to reach a satisfactory final solution coincident with 
the one obtained theoretically (see P] and [10|). 



where ip represent an unknown function (i.e. at the moment 
we have obtained a "partial" solution, in order to reach a 
more satisfactory solution we must take into account the 
Barenblatt criterion) and /3 is: 



1 



3(c.-f 1)(7~1) 

B.2 Calculation of energy density p(t) 

p — p{t, c, Ai^, kj) regarding to the base B, the dimen- 



sional equation of the energy density is: [p] = L ^MT ^ 



A, 



(21) 



B.3 Calculation of radius of Universe f{t). 

f — f{t,c,Ai^,kj) where its dimensional equation is: 

[f]=L 



f (X ct ■ ip 



4(7-l)/3./3 



(22) 



B.4 Calculation of temperature 0{t). 

being its dimensional equation: [6] = 



B. Solutions through D.A. 

We shall calculate through D.A. i.e. by applying Pi- 
Theorem variation of G(t) in function on t, energy density 
p{t), the radius of Universe f{t), the temperature 9{t), the 
entropy S{t), the entropy density s{t) and finally the varia- 
tion of the cosmological "constant" A{t). The dimensional 
method brings us to (see |^ and |^): 

B.l Calculation of G(t) : 

G = G{t, c, A^^ kj) where the dimensional equation of G 
regarding to the base B is: [G] = L'^M'^T'^ . Under this 
circumstances, the application of Pi-Theorem brings us to 
obtain the following dimensionless monomials: 



/ G t 

L 3 

M -1 

\ T -2 1 

TTl 



7 



1 2 + 3tJ 7-1 

11-7 
-1 -2 27-1 y 



ct 



1+/3 



GA^ 



TT2 



(7-l)/3i,/3 



A 



It is observed that the first monomial (tti) is identical 
to the one obtained in the paper (^) for perfect fluids, 
while the second monomial contains information on flow 
viscosity^. These two monomials lead us to the following 
solution: 



G 



ct 



1+0 



A,. 



■ V 



4(7-l)/3./3 



(20) 



■^these remarks, obviusly are valid for all the solutions obtained 
bellow 



A ^ 



,1(1+'^) 



ct^+P 



A 



{-1-^)01,13 



{ct^ - y T-iijj A,^ ^ 

B.5 Calculation of entropy S{t). 
S = S{c,A^a,k^,t) where [S] = 



(23) 



S oc 



{Alaitc)'('- 



4(7-1)/3l/3 



(24) 



B.6 Entropy density s{t). 

s = s{t,c,A^a,k^) where [S] = L-^MT-^O^^ 



(Alay ( ct^+f 
s oc TTTrr. — r ' ^ 



{ctf 



^(7"l)/3^,/3 



(25) 



B.7 Calculation of cosmological "constant" A(t). 
A = A(t, c, At^, k-y) being its dimensional equation [A] 



L- 



1 I ct^+f 



(26) 



C. Different Cases. 



All the following cases that we go on to study now have 
been studied by Arbab (see ^) confirming "j!" his solution 

(0)- 

In obtaining all solutions depending on two monomials 
we shall try to find a solution to the problem expounded by 
means of the Barenblatt criterion (for more details about 
the method used here see S and M). 
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C.l 7 = 1/2 and w = 1/3, Radiation predominance. 

As we pointed out in the introduction the only mod- 
els topologically equivalent to the ones of classic FRW are 
those that follow the law ^ cx ^oP^^^ i-^ 7= 1/2 for its vis- 
cous parameter. In this case we observe a Universe with 
radiation predominance w = 1/3. In order to obtain a com- 
plete solution we shall take into account Barenblatt crite- 
rion since, we have obtained the solutions depending on an 
unknown function ip. In this case the substitution of the 
values of u and 7 leads us to: 



G OC — ; • ip 



1/4.-1/2 



To get rid of the unknown function we apply Baren- 
blatt criterion, for this purpose we need to know the order 
of magnitude of each monomial^: 



10 



-10.59 



7r2 = 



ct 



1/2 



41/4..-1/2 



10 



,2.6 



G(x 




G cx const. 



as we expected in having a model with 7= 1/2. We also 
obtain from this point that fc^ = c^/G. Whit regard to the 
rest of quantities we operate identically finding without 
surprise that: 

fo,t'/^ 0OCi-l/2 Socf 

s OC t^^^ A OC const. 

As we see the model shows the same behavior in the prin- 
cipal quantities as in the classic FRW model with radiation 
predominance. 

Let see, for example, how / has been calculated: Follow- 
ing the same steps as we have seen in the case of calcula- 
tions of G it is observed that: 



focct-^\ 



(27) 



TTl = ^ « 10 

ct 



-2.6 



„^l/2l,l/2 
CI 2.6 
TTS - ~ 10 



, Jct^/^k',/'^ 
J OC ct ' 



A, 



1/4 



-A 



1/2 \ 2 



t^ OC 



/ m : 



GA^ 



t^ 



This solution coincides with the one obtained for a classic 
FRW model with radiation predominance. In any other 
cases fc-y as well as A,^ will have other values to calculate. 

^see the table of numerical values at the end of the text 



C.2 7 = 1/2 and u) = 0. Matter predomineince 

A model with matter predominance w = y topologi- 
cally equivalent to a classic FRW. In this case we find the 
following relationships: 

Regarding to G the solution obtained is (after replacing 
values 7 and co): 



G OC -:— • </? 



A. '^\Ai!%^'') 

as we are working with a model described by matter instead 
of considering energy density we find more convenient to 
consider matter density which becomes a little dimensional 
readjustment in A^ which becomes [AJ\ = M in such a 
way that the solution pointed out above for G is still the 
following law: 



tc-" 



ip 



ctk^ 



1/3 



as in the previous case we apply Barenblatt criterion which 
brings us to: 



TTl = 



gi: 



10 



-1.42 



ctk: 



2 \ 1/3 



7r2 



Au 



10 



-0.47 



TTl = (772)™ /m « -3 

G OC i.e. G OC const. 

In regard to the rest of quantities if we operate as before, 
we get: 



c 1 
poc^ /oc(MG) 3*2/3 A OC const. 

where we have used the equality fc^ = c^/G and we have 
identified A^. with the total mass of Universe M i.e. The 
same behavior has been obtained as in a FRW with matter 
predominance. Let see for instance how we calculate radius 
/: 

For this quantity the obtained solution is: 



, . / ctk 

f QCCt-(p{ 



Barenblatt criterion brings us to: 



2\ 1/3 



f n - / ctk 

TTl = 4 = 10°-" TT2 = 



2\ 1/3 



Ct 



10 



-0.47 



= (7:2)'" /mw-1 /oc(MG)h2/3 
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C.3 7 = 3/4 and = 1/3. An Universe with radiation 
predominance: 

In this case, as /3 = — 1 we find the foUowing solutions: 

ij c { cA^'-v 



as the unknown function does not depend on t we can 
state fearlessly that 



A. 



1/4 



= D = const. 



since c, kj as well as are constant through hypothesis, 
in such a way that 

G cx D't^ 

where D' = Dc^/Ai^. In this case we do not need to re- 
sort to Barenblatt criterion in order to obtain a definitive 
solution. In regard to the other quantities we obtain the 
following behaviors: 



p oz D 



{ctY 



f cx Dct a}I^Q cx D 



A 



1/4 



Ct 



('43 >il/4 

So^D{AlaY'^ so^D ^ " /, A oc ^(ct)-^ 

{cty 

In short, the obtained behaviors are: 

p cx t^'^ f oc t 9 oc t^^ 
S cx const, s cx t^^ A cx t^^ 

this case follows an identical behavior to the one obtained 
in a model described by a perfect fluid with G and A vari- 
ables (see |l|] and Q) showing in this way the generality 
that we have obtained when considering a bulk viscous flu- 
id. 

C.4 7 = 2/3 and uj = An Universe with matter 
predominance. . 

In this case also /3 = — 1, finding the following relation- 
ships as in the previous case: 



tc^ ( ck. 



A 
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that in the previous case leads us to: 

c^t 

Go,D^ 
Au 

where D = ip ^ ^1/4 ^ • Simplifying in the same way, without 
difficulty we reach: 

p (xt~^ / oc t A cx i"^ 

These two last cases are identical to the ones studied in 
references (see Q] and [^). 



V. Not so simple method. 

In this section we will combine dimensional techniques 
with standard techniques of ODEs integration. With the 
dimensional method, we go on to obtain dimensionless 
monomials, which will be replaced in the equations. Thus, 
the number of variables will be reduced in such a way that 
the resulting equation is integrable in a trivial way. We 
study two cases, the first in which we consider div{Tij) = 0, 
while in the other, as we shall see, such hypothesis is not 
needed. 

.5 Considering the condition div{Tij) — 0. . 

In this case we shall pay attention to the equation: 



p' + 3(w -f^ l)pH - <dk^p^H^ + p— 



G' A'c 



'^4 



G 8ttG 







taking into account the relationship div{Tij) = The fol- 
lowing equality is brought up: 



p + 2,{uj + l)pH - 9k-,p''H^=- 



G' A'c' 



'^4 



^ G SttG 



Al 



A2 



The idea is the following: By using D.A. we obtain two 
TT— monomials, which are replaced in the equation, achiev- 
ing a huge simplification of it. On the other hand we in- 
tegrate (^1) and {A2), solving completely in this way the 
problem, this time without Barenblatt. let see. The mono- 

mials obtained are: tti = pk^~~'t-i-^ and 7r2 = Kc^t^ i.e. 



p — akj ^^T-i A = 



where a and d are numerical constants. In a generic way 
the solution is of the following form: p = ak^ "'t^-i if we 



define b = -r-^ then p = ak'it ^ where a 

1— 7 7 



const. G R then 



p' = —hak'^t ^ ^ (paying attention only to the term [Al) 
of the equation) it yields: 

- hak^^t^^-^ + 'i{uj + l)ak^t-''H - 9k^ {ak^^t 

that simplifying it is reduced to: 

9a(''"i) iff - 'dwt-^ff + bt-^f = 

^(^±(u;^-46a-"i)i) 

where w = {lj + 1), if it is defined 
1 



t 



D = 



— j- [w±{w 



Abo?- 



then, the solution has the following form: 
/ = lBt° 





(28) 

(29) 
(30) 

(31) 
(32) 



where I is a certain numerical constant and B is an integra- 
tion constant with dimensions, that can be identified with 
our result by making B = A^k^. 
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Now we shall solve the other term of the equation (the that we replace into the equation that we are treating, re- 



A2). the equation ( 



A'c" 
SttG 



= m 



can be solved suhing: 



in a trivial way if we follow the next results. If we replace 

the monomials tti = phj^"" and tt2 — Kc^t^ in such 

equation the integration of it becomes trivial: 



' G" 







G' = 



6-3 



G{t) 



dc^ 



b~2 



(33) 



where a,d and 5 G R (they are pure numbers). We can 
also observe that this integral needs not be solved since a 
more careful analysis about the number of tt— monomials 
that we can obtain from the equation leads us to obtain a 
solution of the type: 



G — G{k.y^ c, i) 



which brings us to: 



G{t) 



This method, as we have seen, is more elaborated and the 
solution, therefore, finer though coincident with the previ- 
ous one. 

.6 Case in which div{Tij) = is not considered. 

Let see now how we can tackle this problem from 
the D.A. point of view, without imposing the condition 
div{Tij) = 0. The base B as before, is stiU B = {L, A/, T} 
while the fundamental set of quantities and constants this 
time is M = {t, c, fc-y}, with these data we can obtain the 
following monomials from the equation 



p' + 3{uj + \)pH - Qk^p^H^ 



G' A'c 



'^4 



G 8ttG 



(34) 



considering that: 



p = ak^ '' t-i-'^ 



A=-^ (35) 



these two monomials are replaced into the equation, which 
is quite simplified: 



-bak^^t-''-^ 



l)ak^^t-''H 



9k^ {ak^^r^y H"^- 



ak\t-'^ 



dc' 



G AnGt^ 
simplifying this equation, it yields: 



G 



dc^ V 



6-3 



(36) 



(37) 



that along with the field equations (6) and (7) carry us to 
the next set of equations. For example we note that 



-br 



3w 



3c2 




^ 3c2 
that solving it results 



Gr" + — 7 U + 



G' dc^ V 



,2 j-6-3 



3t2 / G Anakb G 



G - gk-'c^t''-^ 



(38) 



where g S M represents a numerical constant. We final- 
ly observe that as in the previous section we could have 
taken into account the three monomials obtained from the 
equation i.e. 



p - afc^t-'' 
replacing them into the equation 

p' + 3{lu + l)pH - 9k,p^H^ + + ^ ^ 

and calculate /, arriving at the same solution obtained in 
the above section i.e. 

/ = IBt^ 

We have proved that it is not necessary to impose the con- 
dition div{Tij) = since it is obtained, in this case, the 
same solution as the one obtain imposing it. 

VI. Conclusions. 

We have studied a cosmological model described by a 
momentum- energy tensor characterized by a fiuid with 
bulk viscosity, in which, furthermore, we have consid- 
ered the constants G and A as functions depending on 
time i.e. as variables and we have imposed the condition 
div{Tij) — 0. We have proved how a suitable use of Dimen- 
sional Analysis enables us to find the solution of such model 
in a "trivial" way. With the "Pretty simple method^ , we 
have obtained two tt— monomials, one of them is the one 
obtained in the case for a perfect fluid ([Q) and the other 
monomial contains the information about viscosity, show- 
ing, in this way, that this type of models is very general 
being able to reproduce the result obtained in the case of a 
perfect fluid. In order to solve the problem we have taken 
into account Barenblatt criterion being able to arrive to ob- 
tain a complete solution of the problem. Standing out that 
our results coincide with the solutions obtained by Arbab 
I. Arbab js). We have shown too that with the "not so 
simple method^ we arrive to solve the problem without ne- 
cessity of impose any condition. We believe, nevertheless, 
that the "simple method^ can be more effective, since, we 
obtain more solutions with it or more complete solutions in 
the sense of flnding in it solutions such as A oc as well 
as A = const, while the "not so simple method^ the only 
solution that is obtained is A oc t"^, but has the drawback 
of depending on Barenblatt criterion i.e. we depend on the 
always insecure numerical data. 



d ^ cH"-^ 



7 





G 


c 


P 


/ 






9 


t 


LO = 1/3 


-10.17 


8.47 


-13.379 


26 


90.62 


13.5 


0.436 


20.25 


= 


-10.17 


8.47 


-26.397 


26 


54 


13.5 
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TABLE I 

The values refer to a logarithmic scale i.e. G ~ lO^-'-''--'-^ 

ETC... MEASSURED IN THE INTERNATIONAL SYSTEM {'m,kg,s}. In 
THE CASE Ul = 0, p REPRESENTS MASS DENSITY WHILE IN THE CASE 
a; = 1/3 REPRESENTS ENERGY DENSITY. 
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